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ABSTRACT 

Main-sequence disruptions of stars by supermassive black holes result in the production of an extended, 
geometrically thin debris stream winding repeatedly around the black hole. In the absence of black hole spin, 
in-plane relativistic precession causes this stream to intersect with itself after a single winding. In this paper we 
show that relativistic precessions arising from black hole spin can induce deflections out of the original orbital 
plane that prevent the stream from self-intersecting even after many windings. This naturally leads to a “dark 
period” in which the flare is not observable for some time, persisting for up to a dozen orbital periods of the most 
bound material, which translates to years for disruptions around black holes with masses ~ lO^M©. When the 
stream eventually self-intersects, the distance from the black hole and the angle at which this collision occurs 
determines the rate of energy dissipation. We find that more-massive black holes (M^ > IO^Mq) tend to have 
more violent stream self-intersections, resulting in prompt accretion. For these tidal disruption events (TDEs), 
the accretion rate onto the black hole should still closely follow the original fallback rate after a fixed delay 
time fdeiay, Afacc(f+fdeiay) = For lower black hole masses (Mh < 10®), we find that flares are typically 

slowed down by about an order of magnitude, resulting in the majority of TDEs being sub-Eddington at peak. 
This also implies that current searches for TDEs are biased towards prompt flares, with slowed flares likely 
having been unidentified. 

Subject headings: black hole physics — galaxies: active — gravitation 


1. INTRODUCTION 

The tidal disruptions of main-sequence stars by supermas¬ 
sive black holes are likely responsible for a few dozen flares 
that have been discovered in the optical (e.g. van Velzen et al. 
2011; Gezari et al. 2012; Chornock et al. 2014; Arcavi et al. 
2014; Holoien et al. 2014; Vinko et al. 2015), UV (Gezari 
et al. 2009), and X-rays (e.g. Komossa et al. 2004; Bloom 
et al. 2011). These tidal disruption events (TDEs) show a 
time-evolution that is consistent with the accretion onto the 
black hole being equal to the fallback rate Mfb set at the 
time of the star’s disruption (Guillochon et al. 2014, hereafter 
GMR14), with the luminosity loosely following a power-law 
decline of (Rees 1988; Lodato et al. 2009; Ramirez-Ruiz 
& Rosswog 2009; Kesden 2012; Guillochon & Ramirez-Ruiz 
2013; Cheng & Bogdanovic 2014). However, it remains un¬ 
clear how these two rates can be so closely related to one an¬ 
other given the difficulty in circularizing material at a distance 
comparable to the star’s periapse distance as dissipation from 
hydrodynamical effects is limited (GMR14). 

General relativistic precession can alleviate this issue if the 
rate of precession is great enough such that the angle of the 
first stream-stream impact is large, and that the intersection 
occurs not too far from the black hole such that material is 
deposited with a semi-major axis comparable to the star’s 
original periapse (Hayasaki et al. 2013; Bonnerot et al. 2015; 
Hayasaki et al. 2015). However, the prototypical disruption 
presented in the literature, a solar-type star being fully dis¬ 
rupted by a IQ^Mq black hole, is decidedly un-relativistic, 
with a periapse distance a factor of ~ 30 times larger 
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than the gravitational radius = GM\i/c^ for a full disrup¬ 
tion (Guillochon & Ramirez-Ruiz 2013). Such a disruption 
does have a stream self-intersection, but at a distance signifi¬ 
cantly larger than the periapse distance, resulting in a viscous 
time fvisc that can be significantly longer than the timescale of 
peak fallback accretion fpeak (Rosswog et al. 2009, GMR14, 
Shiokawa et al. 2015, Bonnerot et al. 2015, Hayasaki et al. 
2015). 

The effect of black hole spin has been considered for the 
effect it has on the initial spread in binding energy of the 
tidal debris (Kesden 2012), and for the effect it has on the 
structure of the debris stream once it has wound around the 
black hole several times (Stone & Loeb 2012; Dai et al. 2013). 
At the moment of disruption, spin can increase (or decrease, 
depending on if the orbit is prograde or retrograde relative 
to the black hole’s spin) the spread in energy by factors of 
a few, leading to flares that potentially peak more (or less) 
quickly than an equivalent Newtonian encounter. In the later 
phases when the stars have been stretched into a thin stream 
(Kochanek 1994), the structure of the stream, and where it 
intersects, can be greatly affected by spin. It is this second 
effect that introduces an additional fixed delay fdeiay between 
the time of the star’s original periapse and when the tidal dis¬ 
ruption is luminous. 

In this paper we focus upon determining fvisc and fdeiay for 
an ensemble of TDEs, and how the values of these parame¬ 
ters both affect the dark period that precedes the start of the 
flare, and its evolution once accretion begins. In Section 2 we 
describe the dynamics of the streams resulting from the tidal 
disruption of a star, and how general relativity (GR) affects 
these dynamics. In Section 3 we describe our Monte Carlo 
approach to constructing a sample of tidal disruption stream 
configurations in order to determine the statistical properties 
of their self-interactions. In Section 4 we present the results 
of our Monte Carlo exercise, and describe the effect of the 
stream dynamics upon fdeiay. Cisc. and the accretion rate Mpeak 
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at time fpeak- We summarize our results in Section 5 and de- cylindrical radius of 
scribe some caveats of our approach, and what our results im- ^ ^ 

ply about the number of identified tidal disruption flares. S = ^ = Wl3rq~^^^, (2) 


2. DYNAMICS OF TIDAL DISRUPTION STREAMS 

Once some material has been liberated from a passing star 
by the black hole’s tides, it is spread in space as a result of 
the range of binding energies it acquired due to the variance 
of the tidal force within the star. This spread is equivalent 
to a spread in velocity of order /3v* relative to the motion of 
the star’s center of mass (Carter & Luminet 1983; Stone et al. 
2013), where v* is the star’s escape velocity, /3 = n/rp the 
star’s impact parameter, is the tidal radius, rp is 

the distance of the star at periapse, and q = Mh/M* is the 
ratio of the black hole’s mass Mh to the star’s mass M». 

Immediately after the encounter, the self-gravity of the 
stream continues to play an important role, confining its width 
to scale as (Kochanek 1994, GMR14), where r is the dis¬ 
tance to the black hole. Self-gravity continues to play a role 
until each part of the stream reaches its apoapse (for the bound 
portion of the stream); at this point the tidal force applied by 
the black hole increases relative to the self-gravity term and 
becomes the dominant component. Because the density of the 
stream continuously drops as the stream stretches in space, 
the parts of the stream that have passed their first apoapse 
are never again confined by self-gravity, and the width of the 
stream scales homologously, i.e. proportional to r. As the in¬ 
ternal energy of the stream is so much smaller than its kinetic 
energy, pressure forces play no role in altering the stream’s 
trajectory, and therefore each part of the stream should closely 
follow a geodesic within the black hole’s spacetime, unless 
the stream intersects with other parts of itself. 

The width of the stream at periapse is determined by the 
ratio of velocity perpendicular to the orbital plane v_l to the 
velocity within the orbital plane vy. The first encounter is 
special; because self-gravity played some role in the stream’s 
evolution, it is somewhat thinner than it would have been oth¬ 
erwise, and thus vx is somewhat smaller than /3v*. Upon re¬ 
turning to periapse a second time, the stream again experi¬ 
ences compression in which approximately 
in specific kinetic energy is converted into internal energy 
through a shock, where /3s is now the effective impact param¬ 
eter for the stream and Cs,s its internal sound speed (GMR14). 
Upon rebound, this energy is converted into kinetic energy in 
the direction perpendicular to the stream’s bulk motion. This 
process repeats every time a part of the stream returns to peri¬ 
apse. 

While simulations and analytical arguments presented in 
GMR14 suggest that the amount of energy dissipated on the 
second encounter is comparable to the first encounter, it re¬ 
mains an open question as to how much energy is dissipated 
on the stream’s subsequent encounters. Simulations that have 
been performed for stars on low-eccentricity orbits (Hayasaki 
et al. 2013; Bonnerot et al. 2015; Hayasaki et al. 2015) suggest 
that the width of the stream does not increase dramatically be¬ 
fore self-intersection occurs (see Figure 13 of Hayasaki et al. 
2013). For this paper we presume that the growth in vx is 
linear with the winding number W, 

y_L=W/3v*, (1) 

where /3v* is the additional spread in velocity imposed each 
time the stream returns to periapse. This leads to a stream 


where W is the number of times the stream has returned to pe¬ 
riapse. Note that this expression is different from that of Dai 
et al. (2013) who assumed that the stream is self-gravitating 
even after repeated windings, however this is likely not the 
case given that the black hole’s gravity becomes dominant 
shortly after the material’s first apoapse approach. Equation 
(1) should be regarded as a lower limit on the stream width; if 
anything, the stream may become wider than what we have 
presumed here given that the shocks formed at the nozzle 
located at periapse may become increasingly violent as the 
stream’s internal sound speed drops relative to its ballistic 
speed (GMR14). As the details of the stream’s width are po¬ 
tentially dependent on cooling and re-ionization (Kochanek 
1994; Hayasaki et al. 2013; Bonnerot et al. 2015; Hayasaki 
et al. 2015), the evolution of the stream’s width remains an 
open question that should be addressed by future study. 

Because the stream cannot occupy a total solid angle that 
is greater than 47r, there is a limit Bjimit on the number of 
windings permitted. 
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where the second expression assumes that Bjimit 3^ 1, al¬ 
though stream self-intersections would likely occur before 
this limit is reached as the relativistic precession principally 
occurs within the orbital plane. As we describe below, the 
first and second windings of the stream almost always inter¬ 
sect in the absence of black hole spin as a result of the apsidal 
precession (Figure 1, left panel), this is significantly smaller 
than the limit presented in Equation (5). 


2.1. Precession About a Non-Spinning Black Hole 

Ignored thus far have been the effects of GR upon the 
stream’s trajectory. Eor non-spinning black holes, the amount 
of precession is dependent upon the ratio of rp to rg, and the 
precession operates prograde to the star’s orbit and entirely 
within its orbital plane, affecting only the argument of peri¬ 
apse U. In Eigure 1 we show the appearance of the stream 
when IT = 5 for non-spinning black holes of four different 
masses, calculated using a freely-available geodesic integra¬ 
tor written for Mathematical (Levin & Perez-Giz 2008). 
Eor the lowest-mass case of = 10^Mq (left panel), pre¬ 
cession is small, whereas Mh > 10 ^Mq shows a significant 
amount of precession, winding more than Itt radians around 
the black hole per periapse passage. 

Dissipation occurs when the stream strikes itself (Rees 
1988; Kim et al. 1999; Kochanek 1994; Ramirez-Ruiz & 
Rosswog 2009; Rosswog et al. 2009; Hayasaki et al. 2013), 
and the amount of dissipation is dependent upon where and 
at what angle the streams intersect, with the specific energy 

3 Available at http : / /goo . gl/FkiUhp. 
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Figure 1. Effect of apsidal precession upon the returning debris stream after the disruption of a star by a non-spinning supermassive black hole. The three panels 
show the structure of the stream after a ^5 = 1 disruption of a 1 Mq star by black holes of four different masses: IO^Mq, IO^Mq, IO^Mq, and with the 

zoom-in focusing upon a region a factor of four times larger than the star’s original periapse. In all panels we show the appearance of the stream a fixed time after 
the disruption, approximately after the most-bound debris has completed ~ 5 orbits, where we have deliberately ignored the effects of stream self-intersection. 
The black disk in all panels is centered about the black hole and has radius rg. 


being dissipated AE^ scaling as 


AP ■ 2 

=-sin 


Q^int 
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(6) 


where r-mt is the distance from the black hole at which the 
intersection occurs, and is the angle at the stream’s self¬ 
intersection point, which can range from zero to tt. Once this 
material has struck itself, it acquires a new binding energy to 
the black hole equal to the amount of energy dissipated (the 
initial binding energy of the debris is negligible). This sets the 
semi-major axis of the resulting disk. 
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which has an orbital period of 



( 8 ) 


If we assume that this disk is subject to the same magneto- 
rotational instability that leads to accretion for steadily- 
accreting black holes, the viscous time (e.g. the time fvisc it 
takes material to accrete) is 


-1 



(9) 


which the fallback rate Mfb(f) would peak, the accretion rate 
onto the black hole as a function of time Macc(f) may or may 
not closely following Mtb(0- In general, Macc(f) = Mfb(t) for 
all times t > fvisc (Shiokawa et al. 2015; Bonnerot et al. 2015; 
Hayasaki et al. 2015); if fvisc ^ tpeak. the flare’s bolometric lu¬ 
minosity should closely follow Mp, at all times, we label such 
events as “prompt.” If fvisc — fpeak. then the accretion onto the 
black hole can be slowed at early times, but the fallback at 
late times and its predicted power-law behaviors (Guillochon 
& Ramirez-Ruiz 2013) may be preserved, an event like this is 
described as being “rise-affected.” In cases where fvisc ^ fpeak. 
the particular shape of Macc is unlikely to survive the circular¬ 
ization process, and the flare is prolonged significantly at the 
expense of reduced peak luminosity, we label such events as 
“slowed.” For slowed flares, the power-law decline is more 
likely to converge to a universal value of f“' (Cannizzo et al. 
1990), and thus valuable information regarding the star’s orig¬ 
inal structure set at the time of disruption may be lost. 

From Figure 1 it is clear that aint for low-mass black holes 
is small, and first occurs near the material’s apoapse, but for 
higher-mass black holes the amount of precession is so ex¬ 
treme that collisions tend to occur closer to periapse and at 
larger angles, with multiple points of self-intersection. The 
only important self-intersection location is the first one, af¬ 
ter this point, the stream is widened significantly and diverted 
from its original path, and the geodesic approximation ceases 
to be appropriate. 


where a is the standard viscous parameter and h is the scale- 
height of the disk (Shakura & Sunyaev 1976). 

As all material in the stream eventually proceeds through 
the self-intersection point, all the mass that eventually ac¬ 
cretes onto the black hole is subject to the same fvisc- Because 
dissipation within the stream is small prior to the stream self¬ 
intersection (GMR14), the specific binding energy of material 
to the black hole dm/de is equal to that set at the time of dis¬ 
ruption. Depending on how fvisc compares to the time fpeak at 


2.2. The Importance of Spin 

When the black hole has spin, the amount of precession 
now depends on the black hole’s spin parameter a and upon 
the inclination of the orbit relative to the black hole’s equato¬ 
rial plane i. Additionally, the precession is no longer confined 
to the orbital plane, resulting in vertical deflection of material. 
This has an important consequence that is not realized around 
non-spinning black holes: Whereas any amount of in-plane 
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of thousands relative to the star’s original size. 

In combination with the “prompt,” “rise-affected,” and 
“slowed” flares described in Section 2.1, the introduction of 
a potentially significant delay caused by the prolonged period 
of time before the stream self-intersects yields six possible 
circularization behaviors, summarized in the bottom six pan¬ 
els of Figure 2. A flare may be for instance “prompt” in the 
sense that its viscous time is short as compared to fpeak, but 
circularization may not begin for years after the star was first 
disrupted. 

Estimating the amount of delay per system is complicated, 
as it depends not only on the physical parameters of the black 
hole and the star (mass, spin, size), but also upon the orien¬ 
tation of the star’s orbit, and on how the stream’s size grows 
with each periapse passage (Equation (2)). A system is likely 
to have a delay if the precession rate is large enough to cause 
the stream to miss itself at least once; this implies that the 
nodal precession per orbit Aw must exceed an angle greater 
than the angular size of the accretion stream upon its return to 
periapse (Stone & Loeb 2012), 


A 

Aw > - : 

rn 


Pq 

27r 


( 10 ) 


We perform a Monte Carlo calculation to determine how often 
this condition is satisfied, and what effects this deflection has 
upon TDE observables. 


Figure 2. Example stream intersections drawn using the Monte Carlo 
method described in Section 3 showing the effects of both apsidal and nodal 
precession. In all panels, the blue ellipses represent parts of the stream that 
do not intersect any other stream, whereas orange ellipses intersect at least 
one other ellipse. Intersection points between different parts of the stream 
are indicated by the small green spheres, whereas the gray spheres show the 
extent of the innermost bound circular orbit (riBCo)- The arrow shows the 
projection of the direction of the black hole’s spin in the image plane, which 
is rotated to remove the precession given to the star upon its first periapse 
passage. The top two panels show a wide (left) and close-up (right) view of a 
stream system in which the most bound material wraps around the black hole 
13 times before intersecting itself. The bottom six panels show prototypical 
examples of the combinations of possible viscous and delay times, with the 
angle and distance of stream self-intersection being denoted in the bottom 
right of each paneP. 

^ A movie showing all stream configurations in the ensemble is available 

at http : / /youtu .be/W3VCdidJ5zg. 

precession would inevitably result in stream intersection in 
the non-spinning case (even for a stream of zero width), the 
nodal precession in a spinning black hole system can deflect 
material out of its original orbital plane such that an intersec¬ 
tion is no longer guaranteed (Dai et al. 2013). 

The main consequence of this is that there can be a signif¬ 
icant delay before stream self-intersection, and thus circular¬ 
ization, can begin. In Eigure 2 we show a small sub-sample 
of the realizations drawn using the method described in Sec¬ 
tion 3 to demonstrate how the nodal precession causes the 
stream to miss itself, with a collision sometimes occurring 
after many windings (Eourteen in the particular example pre¬ 
sented in the top panel). But whereas a viscous time in excess 
of the fallback time can alter the shape of the resulting ac¬ 
cretion function Macc(0. the delay introduced by the stream 
avoiding self-intersection does not affect the shape of 
This is because self-gravity, internal pressure, and dissipation 
at the nozzle-point are all negligible as compared to the bind¬ 
ing energy of the stream to the black hole (Luminet & Marck 
1985), and none of these effects introduce appreciable correc¬ 
tions to Macc(0 even if the stream stretches by factors of tens 


3. MONTE CARLO REALIZATIONS OE TIDAL DISRUPTION 
STREAMS 

The trajectories presented in Eigure 1 were generated in¬ 
tegrating the full geodesic equations around a non-spinning 
black hole, and the code that we utilized is also capable of 
integrating geodesics about spinning black holes. However, 
generation of these curves is computationally expensive given 
that every element of the tidal disruption stream has a differ¬ 
ent binding energy and angular momentum, requiring that all 
points in the stream need to be integrated separately (alter¬ 
native methods that avoid direct integration may be faster, see 
e.g. Dexter & Agol 2009). Additionally, it is difficult to derive 
a robust algorithm for finding when the stream self-intersects; 
the distance between every point on the curve needs to be 
evaluated between every other point, meaning that the com¬ 
putational cost scales with the square of the number of points 
used to resolve the stream. This is further complicated by the 
fact that the stream width scales both with r and W, the deter¬ 
mination of which requires tracking the number of times each 
fluid parcel has returned to periapse. 

As we are interested in generating a sample that’s large 
enough to perform statistics, we want to avoid having to 
solve many differential equations and performing an expen¬ 
sive calculation for each realization. Every fluid element in 
the tidal disruption stream possesses a very large eccentric¬ 
ity ^ as most tidally-disrupted stars originate from 

the sphere of influence where their initial binding energy to 
the black hole is close to zero (Magorrian & Tremaine 1999; 
Wang & Merritt 2004; MacLeod et al. 2012). As a conse¬ 
quence, practically all of the precession resulting from rela¬ 
tivistic effects occurs within a few periapse distances, with 
the trajectory following a Keplerian ellipse beyond these re¬ 
gions (Eigure 1). This means that the general shape of the 
resulting stream is well-approximated by a series of ellipses 
where each ellipse has precessed by some amount relative to 
ellipses corresponding to material with longer periods. 
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Figure 3. Comparison between heuristic prescription (as described in Section 3) and exact solutions for non-spinning (top row) and spinning (bottom row) black 
holes of different masses, where Iq is the stain’s initial inclination relative to the black hole spin plane. Each panel shows two-dimensional projections of the 
exact solution as the dashed black curve, and the heuristic solution with the solid red curve and, with the heuristic curve terminating once a self-intersection has 
occurred, and the exact curve being evaluated such that the number of windings is equal to the heuristic solution. The gray sphere in each panel shows hbco- For 
clarity, the line thickness in this Figure is not equal to the stream’s tme thickness. 


This motivates our method where we represent the tidal 
disruption stream heuristically by a series of unconnected el¬ 
lipses. Such configurations are trivial to generate rapidly, only 
requiring the inclination i and eccentricity e of each ellipse to 
be substituted into the post-Newtonian expressions for preces¬ 
sion (Merritt 2013) 
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( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 


where Auj and A17 are the orbit-averaged nodal and apsidal 
precession terms, and where we have included the de Sitter 
terms (subscript D), Lense-Thirring terms (subscript J), and 
quadrupole terms (subscript Q). 

There are two qualitative behaviors that are lost by this ap¬ 
proach. Firstly, the ellipses do not capture the “zoom-whirl” 
behavior, which manifests when is only a factor of a few 
larger than r^, an example zoom-whirl is visible in the right¬ 
most panel of Figure 1 . However, when spin is introduced, 
the probability of self-intersection within the zoom-whirl re¬ 
gion is small as most of the nodal deflection is occurring there. 
Secondly, because the ellipses are unconnected, we do not ex¬ 
plicitly resolve the small connecting segments between neigh¬ 
boring ellipses. These segments are only ~ long and thus 
represent a small fraction of the stream’s total length, and they 
are somewhat accounted for by the periapse portions of each 


ellipse. Over a large ensemble of realizations, we would ex¬ 
pect that this small difference would not affect the probability 
that a self-intersection occurs for a given set of physical pa¬ 
rameters. 

The system of ellipses is constructed as follows; For each 
orbit of the most-bound material, a new ellipse is created with 
a semi-major axis equal to that of the most-bound material 
fltnb- Ellipses that were previously created are rotated accord¬ 
ing to the precession terms defined in Equations (11) and (12), 
with their semi-major axes being set to amb(lfmb-lY-I- 1)^A^ 
where W^b is the number of windings of the most-bound ma¬ 
terial. 

The orbit-averaged expressions given above are only accu¬ 
rate to 0{v^/c^) and diverge from the true solutions for very 
relativistic encounters. As an example, the expression that is 
expanded via a Taylor series to produce the De Sitter term 
(AH)d becomes imaginary for r < 6rg and is maximally 27r, 
in reality the precession angle can be arbitrarily large (Fevin 
& Grossman 2009). Qualitatively, we And that the precession 
terms slightly underestimate the degree of precession when 
in the weak relativistic regime (leftmost panels of Figure 3). 
Some additional error is also introduced by the orbit-averaged 
rotations being performed in some order. We choose to per¬ 
form the apsidal rotation followed by the nodal rotation, how¬ 
ever the apsidal and nodal rotations are in reality occurring at 
the same time. Both of these simplifications introduce some 
error in the quantities that are expected to be conserved: The 
component of the angular momentum parallel to the black 
hole’s spin axis, and the Carter constant Q (Drasco & Hughes 
2004). We And that the error in L, and (7 is < 10% for 93% 
of realizations, with the greatest deviations occurring for sys¬ 
tems in which r < 6rg. Because the orbital energy for each el¬ 
lipse is set by the number of windings, which depends on the 
time since disruption, orbital energy is perfectly conserved by 
construction. 

Figure 3 shows a comparison between our heuristic ap¬ 
proach and the exact solution for an array of orbits about 
spinning and non-spinning black holes. From this figure the 
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heuristic solution can be seen to be fairly close to the exact 
solution forMh < with significant deviations begin¬ 

ning to arise for more massive black holes. However, de¬ 
spite the differences in orbital shape, the location and angle 
of the first stream self-intersection is very similar for the ex¬ 
act and heuristic approaches, with the exact solutions tending 
to produce self-intersections that occur slightly closer to the 
black hole and at slightly larger angles. The biggest qualita¬ 
tive disagreement comes for very massive, non-spinning black 
holes (upper right panel of Figure 3) in which the orbit under¬ 
goes a zoom-whirl behavior, resulting in the stream striking 
itself very close to the event horizon. However, as we find 
below, most black holes in which such extreme in-plane pre¬ 
cession exists are likely to be spinning fast enough to deflect 
the stream out of its original orbital plane. 

With the creation of each ellipse, a root-find is performed 
between all pairs of ellipses (labeled i and j) to And the mini¬ 
mum of the following expression, 

X>(xi,xj)= |xi-xj|-5(|x,|)-5(|xj|), (18) 

where S is calculated using Equation (2). If I? < 0, the two 
ellipses are considered to be overlapping and a stream self¬ 
intersection has occurred. In the event that two ellipses in¬ 
tersect at multiple locations or that more than two ellipses in¬ 
tersect (e.g. the 2‘’‘^ winding intersects with both the 4* and 
5* windings), we only consider the intersection that occurs 
chronologically first, as presumably the stream would be de¬ 
flected significantly at the location of this first intersection. 
The viscous time fyisc is then calculated using Equation 9, with 
Tint and Oint being determined at the point where the two el¬ 
lipses overlap. The delay time is calculated by the number of 
windings prior to the first self-intersection, 

fdelay = WP„b, (19) 

where Pmb is the orbital period of the most-bound material. 

We draw random parameters for 4,096 systems assum¬ 
ing a flat distribution in logjgMh with 5 < logjgMh < 8, 
a Kroupa distribution for M* (Kroupa et al. 1993) with 
-1 < logjgM* < 2, a flat distribution in black hole spin 
a, an impact parameter distribution that scales as /3“^ (Rees 
1988), and random initial inclinations i relative to the black 
hole’s spin plane. Stellar radii are calculated using the zero- 
metallicity main-sequence mass-radius fitting formula of Tout 
et al. (1996). We check that the periapse of each star lies out¬ 
side the black hole’s IBCO (Bardeen et al. 1972), if not we 
redraw (3 and i, keeping other parameters fixed, until the peri¬ 
apse lies outside this region. 

4. DELAYS AND SLOWDOWNS FOR TIDAL DISRUPTION EVENTS 

Our Monte Carlo exercise allows us to calculate the dura¬ 
tion of the dark period fjeiay and the viscous time fyisc for a 
random ensemble of TDEs. We And in general that flares 
with a significant f^eiay tend to have their first stream self¬ 
intersection once the most-bound material has wrapped sev¬ 
eral times around the black hole, with a half-dozen wind¬ 
ings being common, and seventeen windings being the maxi¬ 
mum found in our ensemble of realizations (Eigure 4). Usu¬ 
ally the stream self-intersection occurs between neighboring 
windings, e.g. the 5* winding is likely to have an intersection 
with the 6* winding, but unlikely to with the 12* winding. 
This implies that that ratio of densities between the two wind¬ 
ings are likely to be similar to one another, and thus these self¬ 
intersections should lead to significant dissipation. Events 
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Figure 4. Density histogram of realizations as a function of the winding 
numbers of the two parts of the stream that intersect one another, where Wi 
and W 2 refer to the winding number of each ellipse. As a single ellipse cannot 
intersect itself, Wi = W 2 is not permitted. 
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Figure 5. Dependence of Lisc/lpeak on the amount of apsidal (AD) preces¬ 
sion per orbit. The background coloring corresponds to the circularization 
behavior described in Section 4, with the color of the points corresponding to 
logjQMh as indicated by the colorbar. 

with the smallest fjeiay values tend to be those in which the 
U' and 2"“* windings intersect, in other words these events be¬ 
gin circularizing once the most bound material completes a 
single orbit around the black hole. 

In Eigure 5 we show how fvisc and fjeiay relates to the time 
of peak accretion fpeak. where we use fitting relations (Guillo¬ 
chon & Ramirez-Ruiz 2013, 2015) to determine fpeak- For the 
purposes of calculating fvisc, we need to assume the param¬ 
eters that govern the viscous evolution of the accretion disk 
once it has formed. Eor the figures and fitting formulae in 
this paper we select values appropriate to a thick disk, a = 0.1 
and (h/r)^ = 0.1, yielding a viscous time that is 100 times 
longer than the orbital period (this assumption is relaxed in 
Section 5). We consider an event to be “prompt” (green color 
in the figures) if 3fvisc < fpeak, to be “rise-affected” (yellow) if 
fvisc < fpeak < 3fvisc, and to be “slowed” if fpeak < fvisc- 

Eigure 5 shows that Uisc/fpeak is closely related to the 
amount of apsidal precession per orbit, and because this gen¬ 
erally increases with increasing Mh as r^/rg oc ' , the vis¬ 
cous times tend to be shorter for more-massive black holes. 
There are exceptions even for large values of AH which still 
have large fvisc values, these are frequently systems with long 
delay times that happen to have stream self-intersections with 
small intersection angles aint (e.g. lower right panel in Fig- 
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Figure 6. Dependence of ^deiay upon the black hole’s spin parameter a. Each 
point shows an individual realization, with the color of the points correspond¬ 
ing to Mh as indicated by the color bar. The colored diamonds show the 
mean values within a range of spin parameter indicated by their eiTor bars, 
where we have combined the prompt and rise-affected categories into the yel- 
low/green diamonds, and the slowed events are shown by the red diamonds. 


ure 2). Because some stream self-intersections occur at nearly 
zero Oint resulting in little dissipation, fyisc is set to be the mini¬ 
mum between Equation (9) and the time required for the com¬ 
bined cross section of all parts of the stream to cover the full 
47r steradians, 

^cover “ EjYibfklimiti (20) 

where VTiimit is calculated using Equation (5) and Pmb is the 
period of the most-bound material. This typically caps fvisc to 
be no more than a factor of 100 times greater than fpeak- 

The time of delay fdeiay between when the disruption of the 
star occurs and the accretion of matter onto the black hole 
is shown in Eigure 6 as a function of the spin parameter a. 
While there is a lot of scatter at a given a, as one might ex¬ 
pect given the range of initial inclinations relative to the spin 
plane, there is a clear reduction in mean values of fdeiay for 
events with a < 0.2. This suggests that there is a minimum 
a required to deflect the stream enough to avoid striking it¬ 
self in the first few windings, but once above this threshold, 
the spin of the black hole does not have a strong influence on 
fdeiay Because the observed spins of supermassive black holes 
tend to be greater than 0.5 (Reynolds 2014) as opposed to our 
assumption of a flat distribution of a, this suggests that most 
black holes possess enough spin to induce the deflections nec¬ 
essary to cause a delay. As a < 0.2 black holes constitute 
only 20% of our sample, we And that their exclusion does not 
alter the statistical properties of the ensemble, and thus our 
results should be directly applicable to the observed, higher- 
spin sample. We also And that fvisc is very slightly shorter for 
larger a, suggesting that the main effect of a is to increase the 
amount of time between disruption and circularization, but 
has little effect on the circularization duration. 

In Eigure 7 we show a scatter plot of fdeiay as a function of 
Mh, where the color-coding corresponding to its circulariza¬ 
tion behavior. With our assumptions on the viscous time, we 
And that there is a change in the fraction of prompt events at 
a black hole mass of with the majority of events be¬ 

ing slowed below this mass and the majority being prompt or 
rise-affected above it. If the viscous time were to be shorter 
than what we have assumed, the transition mass would move 
to lower black holes masses and more events will be consid¬ 
ered prompt, whereas a longer viscous time will slow down 
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Figure 7. Length of the delay period /delay as a function of black hole mass 
M^. The points are color-coded according to how the viscous time /vise com¬ 
pares to /peak’ the time of peak for the fallback rate. Prompt events are rep¬ 
resented by the filled green points, events with an affected rise by the filled 
yellow points, and slowed events by the open red points, these classes are 
defined in Section 4. The histograms on the right and top sides of the figure 
show the fraction of events belonging to each category. The solid lines show 
a least-squares fit to only the prompt flares, whereas the dashed line shows a 
fit to the slowed events only. 

prompt events occurring around black holes of higher mass. 

We perform a least-squares fits to fdeiay to determine the 
trends in fdeiay as a function of M^. Eor events that are slowed, 
we And that the delay time scales sub-linearly with the black 
hole’s mass. 


tdelay.slow “ 0.92 


/ Mh 


VIO^Mg 


0.78 


yr, 


( 21 ) 


whereas we And that for events that are prompt follow a bro¬ 
ken power-law behavior, with events falling below the cutoff 
black hole mass having fdeiay scaling super-linearly with Mh, 


tdelay,fast ' 


I 0.93 ( 

^ Mk ) 

, lO’Mo ) 

] 0.97( 

' Mk \ 

, lO’Mg ) 


1.44 


yr, Mh < 1O^ °M0 
0.82 ■ (22) 
yr, Mh > 10^ °Mo 


Because no accretion takes place prior to circularization, 
this delay period is dark relative to the peak flare luminos¬ 
ity, with the only light likely coming from recombination of 
the stream as it expands, with luminosity L < 10"^' ergs s“' 
(Kasen & Ramirez-Ruiz 2010). This delay has implications 
for the time one expects to detect signatures of a disruption 
arising from the moment of maximum compression when the 
star first passes periapse, such as shock break-out (Brassart & 
Luminet 2008; Guillochon et al. 2009) and gravitational wave 
emission (Kobayashi et al. 2004; Rosswog et al. 2009; Stone 
et al. 2013). Eor TDEs without a significant delay, these sig¬ 
natures would be expected to be detectable weeks to a few 
months prior to the peak of the accompanying accretion flare. 
Delayed events, on the other hand, can show these compres¬ 
sion signatures years prior to the flare. 

Eigure 7 also shows that slowed and prompt events co¬ 
exist for black holes of the same mass, indeed there are a 
small fraction of events that are prompt even forMh = lO^M©. 
Around low-mass black holes (Mh < lO^M©) the tidal radius 
is significantly larger than for main-sequence star disrup¬ 
tions. Therefore, prompt flares about these black holes require 
deeply plunging encounters, with /3 > 5 being typical. 
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Figure 8. Histograms showing the ratio of the peak to Eddington accretion rates (left panel) and (right panel) for the initial material returning to 

pericenter Mn, (dashed black line) and the accretion rate onto the black hole Macc (solid black line) for the three classes of TDEs defined in the text; color-coding 
matches that of Figures 6 and 7. The three sub-panels show three different ranges of . In the left panel, the dotted vertical line shows where Mpeak = ^Edd ■ 


Because the Eddington limit LEdd = 47rGMhmpc/crT for a 
black hole scales with Mh, TDEs about less-massive black 
holes tend to have peak fallback rates that are a few orders 
of magnitude greater than this limit, whereas more-massive 
black holes are usually sub-Eddington near peak. The slow¬ 
ing of flares around lower-mass black holes reduces their peak 
accretion rate relative to the fallback rate, and prolongs the 
flares such that they peak at a later time, whereas a large frac¬ 
tion of disruptions about higher-mass black holes are prompt, 
meaning that their accretion rates are not strongly affected by 
circularization. This has the net effect of reducing the peak ac¬ 
cretion rate such that super-Eddington events are significantly 
more rare than if all TDEs were prompt. 

In Eigure 8 we show Mpeak/A^Edd and fpeak, where we have 
presumed that the accretion efficiency /i = 1 - -^1 -2/3risco 
and risco is calculated as in Bardeen et al. (1972), and we use 
the fitting relations of Guillochon & Ramirez-Ruiz (2013) to 
determine Mpeak and fpeak- By assuming that the peak accre¬ 
tion rate is reduced by a factor that is inversely proportional 
to fvisc/tpeak^ we find 56% of TDEs are super-Eddington once 
this slowdown effect is included for Mq < Mh < 10 ®Mq 
versus 82% if all TDEs were prompt, 34% for 10® Mq < 
Mh < IQPMq (versus 49%), and 2.6% for 10 ^Mq < Mh < 
10*^ Mq (versus 4.0%). Because more-massive black holes 


have higher Eddington luminosities, sub-Eddington events 
from higher-mass black holes will typically be more luminous 
than flares from lower-mass black holes. As the typical Ed¬ 
dington ratio only scales weakly with Mh for Mh <10^ (as 
shown in the top two left panels of Eigure 8), and the Edding¬ 
ton accretion rate scales as Mh, TDEs from black holes from 
1O®-1O^M0 end up with similar peak accretion rates at peak, 
with a median value of Mpeak = 10“' ^MQyr“* over this mass 
range. Assuming a radiative efficiency of 10%, this translates 
to a peak luminosity of 10“'"' ergs s“', comparable to the in¬ 
ferred bolometric output of most observed TDEs. Because 
photons are largely trapped once the accretion rate exceeds 
Eddington (Shakura & Sunyaev 1973), the total disk luminos¬ 
ity is potentially limited to the Eddington luminosity (though 
see Jiang et al. 2014), which would suggest that detections in 
those bands are further biased towards sub-Eddington events 
(De Colle et al. 2012). 

This slowdown also affects the timescale over which the 
flare evolves, which is shown in the right panel of Eigure 8. 
Eor lower-mass black holes, the slowdown is quite dramatic 
in most cases, causing some events to peak over timescales 
of several years. Given that most of the claimed opticaEUV 
flares that were observed during maximum light peaked over a 
period of weeks or months (e.g. Gezari et al. 2012; Chornock 
et al. 2014; Arcavi et al. 2014; Holoien et al. 2014; Vinko 
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et al. 2015), and given that no tidal disruption candidates have 
been identified peaking on decade-long timescales, it is very 
possible that these slowed flares are mistaken for longer-term 
AGN activity. This may be the reason why tidal disruption 
surveys, which were originally predicted to detect dozens of 
flares per year (van Velzen et al. 2011), have been discovering 
flares at a rate approximately one tenth that predicted from 
disruption rate estimates (Wang & Merritt 2004; van Velzen 
& Farrar 2014; Stone & Metzger 2014). Moving to higher- 
mass black holes, the slowdowns affect a smaller fraction of 
the total, with the prompt fraction peaking on timescales of a 
few months to a year. 

Because the prompt flares peak on shorter timescales, and 
because their bolometric light curves are likely to follow the 
original fallback rate Mfb (thus preserving the behavior), 
it is this sub-population of the TDEs that are likely to have 
been successfully identified as TDEs. The remainder, which 
are slowed significantly and are likely to have a time depen¬ 
dence that differs from the fallback rate predictions, may have 
been missed by current surveys. This argument also applies to 
the jetted TDEs (e.g. Bloom et al. 2011; Zauderer et al. 201 1), 
which require super-Eddington accretion rates to power the 
jet (De Colle et al. 2012) and rapid accretion timescales to 
explain the observed decay rate (Tchekhovskoy et al. 2014; 
Kelley et al. 2014). 

5. SUMMARY 

In this paper, we performed Monte Carlo realizations of 
tidal disruption streams to determine their structure prior to 
the onset of circularization. We have shown that the effects of 
GR are important when considering the evolution of the de¬ 
bris stream resulting from the tidal disruption of a star by a 
supermassive black hole. When accounted for, we find that 
GR has the following effects: 

• When GR effects are weak and the precession angles 
are small, the stream typically intersects itself far away 
from periapse. This results in a long viscous time that 
can potentially affect the accretion rate Mace and the 
time of peak accretion fpeak (GMR14), a result that has 
been recently verified by numerical simulations (Sh- 
iokawa et al. 2015; Bonnerot et al. 2015; Hayasaki et al. 
2015). The slowdown can result in fpeak increasing by 
a factor of ^100, which also reduces Mace by a simi¬ 
lar factor relative to the fallback rate Mfb- This mostly 
affects disruptions around black holes with masses less 
than a few lO^M©. 

• For more-massive black holes, GR effects are signif¬ 
icantly stronger, resulting in stream self-intersections 
that occur at large angles and closer to the star’s original 
periapse. If the black hole spin is small (a < 0.2), the 
accretion rate onto the black hole typically follows the 
fallback rate, with no significant delay between the re¬ 
turn of the most-bound material and accretion onto the 
black hole. 

• If a black hole is both massive and possesses at least 
a mild spin {a > 0.2), the stream self-intersection oc¬ 
curs after the most-bound material has wound around 
periapse multiple times. Because very little energy 
is dissipated in the stream in the absence of stream 
self-intersections, the specific binding energy of the 


material to the black hole is unchanged until circu¬ 
larization begins, at which point the accretion rate 
will equal the fallback rate with a fixed time delay, 

4^acc(^ “t" ^delay) “ 4/fb(t). 

• Because tidal disruptions around less-massive black 
holes tend to be slowed down, and because more- 
massive black holes have higher Eddington luminosi¬ 
ties, the typical Eddington ratio for a tidal disruption 
at peak luminosity is less than one, with most super- 
Eddington flares originating from black holes with 
masses < IO’Mq. 

• TDEs that are significantly slowed by long viscous 
times likely peak on timescales of several years, with 
bolometric luminosities decaying at a rate more shal¬ 
low than If these slowed events were observed in 
the past, these differences might have prevented their 
identification as TDEs, as current surveys have focused 
on more-rapidly peaking transients. This may explain 
the relative dearth of events as compared to theoretical 
expectations. 

The advantage of our approach is that we can generate a 
large ensemble of tidal disruption streams at little computa¬ 
tional expense, but is heuristic in the sense that we made ap¬ 
proximations to the precession terms using an orbit-averaged 
post-Newtonian formalism that introduces some error, as 
shown in Figure 3. A Monte Carlo approach that utilizes 
the exact solutions is computationally expensive for three rea¬ 
sons: First, we have found that the exact solutions take a few 
minutes to compute on a single CPU per realization, as com¬ 
pared to the near-instantaneous computation time (a fraction 
of a second) of the heuristic used here. Second, the prob¬ 
lem of finding when and where a time-dependent parametric 
trajectory comes within a certain distance of itself is a three- 
dimensional root-finding problem, with the two positions on 
the curve and the time since disruption as free parameters. 
In our heuristic approach, time is discretized by the number 
of windings, so the self-intersection calculation is reduced 
to finding the solutions of Wi(Wi - 1) two-dimensional root 
solves. Third, an even distribution of particles in energy does 
not result in a well-sampled elliptical spiral as the particles 
spend the majority of their time near apoapse. An accurate 
representation of the curve requires a careful selection of par¬ 
ticle energies to represent the debris, with the optimal sample 
being different for every time t. The combination of these 
three issues mean that for every f the full geodesics must be 
integrated from time zero using either a unique optimal sam¬ 
pling of particles (which takes minutes), or a large number 
of evenly spaced particles that guarantees the elliptical spiral 
will be well-sampled at all times. In either case, the total com¬ 
putation time is hours per realization, and is not suitable for 
generating an ensemble of thousands of systems. It is likely 
that alternative geodesic determinations that are faster to com¬ 
pute (such as Dexter & Agol 2009) are required to make the 
exact Monte Carlo feasible. 

We have asserted that the width of the debris stream grows 
via additional kinetic energy deposited each time the stream 
returns to periapse, but the evolution of the stream width as 
a function of time may be strongly influenced by cooling, re¬ 
combination, shock heating, and magnetic fields. Cooling and 
recombination may cause the stream to fragment, resulting 
in the formation of discrete clumps within the stream (Roos 
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Figure 9. Fraction of events that are prompt, rise-affected, and slowed as a 
function of (h/ry'^ below (top panel) and above (bottom panel) the break 
in delay times at IO^Mq. The vertical dashed line shows the fiducial value 
of oT^Qilry^ = 100 used to generate the previous plots and fitting relations 
presented in this paper. 


1992; Guillochon et al. 2014), but will not affect the vertical 
spread in velocity as the stream will not be self-gravitating af¬ 
ter its first return to periapse. However, heating of the stream 
through an oblique shock at periapse (Matzner et al. 2013; 
Salbi et al. 2014) may increase this vertical spread, this would 
increase the typical stream widths and would result in streams 
that self-intersect after fewer windings than what is presented 
here. Lastly, the magnetic fields that permeate the star likely 
thread through the stellar debris and potentially become dy¬ 
namically relevant once the gas pressure drops sufficiently, 
this additional pressure support may also increase the width 
of the stream. 

We have assumed that material circularizes when two 
streams intersect, but this may not be the case if the density 
of the two streams are wildly different from one another, as a 
denser portion of stream may uneventfully pass through a por¬ 
tion of signihcantly lower density. But because neighboring 
windings tend to be the ones we find intersecting (Figure 4), 
the density ratio between the two colliding streams is expected 
to be similar. One improvement that could be made to our cal¬ 
culation would be to map dm/de to each ellipse to determine 
their relative densities. 

Our results are also contingent upon the viscous evo¬ 
lution of the disk, and in this paper we assumed that 
a“'(/!/r)“^ = 100, a value typical of thick disks about 

steadily-accreting AGN. However, viscous times that are both 
shorter and longer than this hducial value are possible. In 
thin accretion disks, a ^ 0.1 and (h/r)^ ^ 10“^ resulting in 
a~^{h/ry^ = 10^ (King et al. 2007), whereas ADAF flows can 
have a = 0.3 (Narayan & McClintock 2008) and (h/r)^ ~ 0.25 
(Yuan & Narayan 2014) resulting in a~^(h/r)~^ = 10. In Fig¬ 
ure 9 we vary the value of and calculate the frac¬ 


tion of events that are prompt, rise-affected, and slowed. For 
disks that are more viscous than our fiducial assumption, we 
find the number of prompt flares increases with decreasing 
a“'(/r/r)“^, although this ratio is maximally 55% (95%) be¬ 
low (above) the break in delay time at 10^Mq even for a 
highly viscous disk where a~'^(hjr)~^ = 1. For less viscous 
disks, the prompt fraction drops to near zero for black holes 
above the delay time break, and accounts for only 5% of 
flares below the break. While it is not clear if the value of 
a~^{h/r)~^ depends on the fallback rate exceeding Eddington 
as it does for steadily-accreting AGN, the strong dependence 
of the prompt fraction on this value (especially for massive 
black holes) suggests that flares require < 100 to 

be observable. 

Despite these uncertainties, our work presented here 
demonstrates that black hole spin plays an important role in 
the circularization and eventual accretion of matter by su- 
permassive black holes that have recently tidally disrupted a 
main-sequence star. These effects may also be important for 
other stellar disruptions in which is comparable to r^, such 
as the disruption of a white dwarf by an intermediate mass 
black hole (Rosswog et al. 2009; Krolik & Piran 2011; Haas 
et al. 2012; MacLeod et al. 2014), or the disruption of a gi¬ 
ant star by a very massive (> IQ^Mq) black hole (MacLeod 
et al. 2012, 2013). Future (magneto-)hydrodynamical simu¬ 
lations should aim to establish how the width of the stream 
evolves through the repeated periapse encounters made pos¬ 
sible by the deflection induced by black hole spin, and the 
viscous evolution of the debris once it has circularized. 
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